Optimization of target biochar for the adsorption of target heavy metal ion

The purpose of this work is to study the pyrolysis conditions of target biochar suitable for target heavy metal ion, to characterize the optimized target biochar, and to study the adsorption performance of biochar. With Cu2+ and Zn2+ as the target pollutants, the pyrolysis conditions involved in the preparation process as pyrolysis temperature, pyrolysis time, and heating rate were evaluated and optimized from Box–Behnken Design (BBD), response surface methodology (RSM) and desirability function, the optimized pyrolysis conditions of target biochar for Cu2+ (Cu-BC) and Zn2+ (Zn-BC) were obtained. The optimum pyrolysis parameters for Cu-BC and Zn-BC were pyrolysis time of 3.09 and 2.19 h, pyrolysis temperature of 425.27 and 421.97 °C, and heating rate of 19.65 and 15.88 °C/min. The pseudo-second-order kinetic and Langmuir isotherm model proved to be the best fit for the equilibrium data, with a maximum adsorption capacity (Qmax) fitted by Langmuir model were 210.56 mg/g for Cu2+ by Cu-BC and 223.32 mg/g for Zn2+ by Zn-BC, which were both higher than the Qmax of unoptimized biochar (BC) for Cu2+ (177.66 mg/g) and Zn2+ (146.14 mg/g). The physical properties, chemical structure, surface chemistry properties of Cu-BC and Zn-BC were characterized by Zeta potential meter, Scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDX), Fourier-transform infrared spectroscopy (FTIR), and X-ray diffraction (XRD). This study puts forward a new perspective for optimizing target biochar production for special environmental application.


Preparation of biochar.
Cut the airing water hyacinth into 1-2 cm segments, and dried them in drying oven at 70 °C until constant weight. The dried water hyacinth was put into the high temperature resistant quartz boat, and then placed in a tubular furnace. The preparation method is shown in the previous study 17,18 . The pyrolysis parameters (pyrolysis temperature, pyrolysis time and heating rate) were set according to the experimental requirements, and the whole process was carried out under the protection of nitrogen. The specific experimental operations of the three preparation parameters are as follows: 1) Pyrolysis temperature. The pyrolysis time and heating rate were fixed at 2 h and 20 °C/min, respectively. The pyrolysis temperature was selected as 200, 300, 400, 500, 600 and 700 °C, respectively.
The pyrolysis temperature and heating rate were fixed at 400 °C and 20 °C/min, respectively. The pyrolysis time was selected as 1, 2, 3, 4 and 5 h, respectively.
3) Heating rate. The pyrolysis temperature and pyrolysis time were fixed at 400 °C and 2 h, respectively. The heating rate was selected as 5, 10, 15, 20, 25 and 30 °C/min, respectively.

Adsorption experiments of target biochars.
(1) Screening of central value of preparation conditions for target biochar.
In order to investigate the central value of preparation conditions on the adsorption of target pollutants by biochar, 2.0 g/L of biochars were reacted with 20 mg/L of Cu 2+ or Zn 2+ solutions for 240 min on the oscillator SHA-CA at 150 rpm/min at 25 °C. The Cu 2+ and Zn 2+ stock solutions were prepared with Cu(NO 3 ) 2 ·3H 2 O and Zn(NO 3 ) 2 ·6H 2 O, respectively.
With Cu 2+ or Zn 2+ adsorption capacity of biochar as response, the preparation conditions of the target biochar were optimized by using the Box-Behnken Design (BBD), and the interaction between various conditions were explored by response surface methodology (RSM). As the operation conditions, the pyrolysis time (X 1 ), pyrolysis temperature (X 2 ), and heating time (X 3 ) were investigated by RSM, and the Cu 2+ or Zn 2+ adsorption capacity of biochars was selected as the response variable (Y Cu or Y Zn ). The experimental results were analyzed and fitted to a quadratice quation by using the Design Expert 10.0 software. In the optimization experiment, the concentration of heavy metal ions changed to 50 mg/L, and other experimental conditions were same as the preliminary test. The biochars prepared under optimized preparation conditions were called Cu-BC (for Cu 2+ ) and Zn-BC (for Zn 2+ ).
(3) Adsorption kinetics of target biochar for target heavy metal ions.
Pseudo-first-order model: Pseudo-second-order model: q e and q t are adsorption capacity of Cu 2+ or Zn 2+ adsorbed by biochar at equilibrium and t time, respectively (mg/g); k 1 : the pseudo-first-order rate constant (1/min); k 2 : the pseudo-second-order rate constant (g/mg min). (4) Adsorption isotherm experiments of target biochar for target heavy metal ions.
Adsorption isotherm experiments for target heavy metal ions were conducted by dosing 0.2 g of target biochars in 100 mL of solutions with a range of Cu 2+ and Zn 2+ concentrations (10,20,50,100,200, 500, 800 and 1000 mg/L).
The pH values of the heavy metal solutions in all batch experiments were adjusted to about 5.5 by 0.1 mol/L HNO 3 and NaOH solutions. The concentrations of Cu 2+ and Zn 2+ after adsorption were determined by Shimazu ICPE-9000 inductively coupled plasma emission spectrometer. All experiments were repeated in triplicate. The equilibrium data were analyzed by the Langmuir and Freundlich models, and the two models were shown in Eqs. (4) and (5) 17,18 .
Langmuir model: Freundlich model: q e or Q e : adsorption capacity of Cu 2+ or Zn 2+ by biochar at equilibrium (mg/g); C e : concentration of Cu 2+ or Zn 2+ at equilibrium (mg/L); Q m : the Langmuir maximum adsorption capacity (mg/g); K a : Langmuir constant (L/mg); K F : the Freundlich constant; 1/n: adsorption strength of the system.
Characterization of target biochars. The characterization methods of biochar were detailed in previous research 17,18 . Scanning electron microscope microscopy and energy dispersive spectroscopy (SEM-EDX) analysis was conducted using S-4800 scanning microscope and energy dispersive X-ray spectroscopy (Hitachi, Japan) to examine the surface topography and the element types and contents of biochar. Specific surface area of biochar was determined using a NOVA 2000e analyzer. Nitrogen adsorption-desorption was carried out in 77 K liquid nitrogen. The specific surface area and pore size distribution of the adsorbed material were calculated according to the peak area of adsorption and desorption. The functional groups of biochars were analyzed by an SHIMADZU IRPrestige-21 transform infrared spectrometer (Shimadzu, Japan). 1 mg sample powder was mixed with KBr at the mass of 1:200 and ground evenly. After pressing, the samples were put into the sample chamber and scanned in the wavelength range of 400-4000/cm. The crystallinity changes of the three materials were characterized by X-ray diffraction (XRD) with Bruck-D8 series X-ray (powder) diffractometer (Bruck, German). 0.2 g of sample powder was taken and put into the sample chamber after tablet pressing, and scanned in the range of 10°-80°. (2) q t = q e 1 − exp(−k 1 t)  Fig. 1a-c. Figure 1a showed the pyrolysis time effects on the adsorption capacity of Cu 2+ and Zn 2+ by biochars. It can be seen that the adsorption capacity of Cu 2+ and Zn 2+ by biochars increased first and then decreased with the pyrolysis time prolonging. According to the Fig. 1a, the optimal pyrolysis time of biochar corresponding to Cuadsorbed biochar is 3 h, and the maximum adsorption capacity of Cu 2+ is 9.87 mg/g, while the optimal pyrolysis time of biochar corresponding to Zn-adsorbed biochar is 2 h, and the maximum adsorption capacity of Zn 2+ is 9.79 mg/g. According to the above experimental results, the pyrolysis time of 3 h and 2 h were selected as the center values for Cu-adsorbed and Zn-adsorbed biochar, respectively. So, in the optimization experiments, the low, medium and high levels for Cu-adsorbed biochar were 2 h (− 1), 3 h (0) and 4 h (+ 1), and for Zn-adsorbed biochar were 1 h (− 1), 2 h (0) and 3 h (+ 1). The codes − 1, 0, and + 1 indicate low, medium and high levels of preparation parameters.
The effect of pyrolysis temperature on the adsorption property of biochars is displayed in Fig. 1b. From the Fig. 1b, the adsorption capacity of biochars on Cu 2+ and Zn 2+ improved with the increase of pyrolysis temperature. When the pyrolysis temperature increased to 400 °C, the adsorption capacity of biochar for Cu 2+ and Zn 2+ reached 9.68 mg/g and 9.79 mg/g, respectively. Continued to increase the pyrolysis temperature to 500 °C, the adsorption capacity decreased slightly, while the pyrolysis temperature is further increased to 600 °C and 700 °C, the adsorption capacity of Cu 2+ and Zn 2+ were higher than that of 500 °C, but lower than that of 400 °C. Therefore, 400 °C was chosen as the optimal temperature for the production of biochar from water hyacinth adsorbing Cu 2+ and Zn 2+ , and the three levels of pyrolysis temperature in BBD were 300 °C (− 1), 400 °C (0) and 500 °C (+ 1). Figure 1c represented the effect of heating rate on the adsorption capacity of biochars for Cu 2+ and Zn 2+ . With the increase of heating rate, the lignocellulosic structure of water hyacinth gradually decomposed, and the volatiles precipitated promoted the formation of biochar pores. When the heating rate continues to increase, the volatiles escape sharply, the residual pores become larger or some products block the pores, as a result, the specific surface area of the prepared biochar decreased, and the adsorption performance of the biochar is reduced. From the Fig. 1c, when the heating rate was 20 °C/min, the adsorption capacity of biochar on Cu 2+ was the highest, reaching 9.83 mg/g, while for Zn 2+ , the optimal heating rate was 15 °C/min, the adsorption capacity of Zn 2+ was 9.84 mg/g. Therefore, 20 °C/min and 15 °C/min were selected as the center values for preparing biochar to adsorb Cu 2+ and Zn 2+ in the BBD, respectively. The three levels of heating rate in the BBD were 15 °C/min (− 1), 20 °C/min (0) and 30 °C/min (+ 1) for Cu-adsorbed biochar, and 10 °C/min (− 1), 15 °C/min (0) and 20 °C/min (+ 1) for Zn-adsorbed biochar.
Optimal preparation conditions of target biochars. Table 1 summarized the actual and coded values of the preparation parameters for biochars. Meanwhile, a total of 17 runs of randomized BBD experiments were designed, and the corresponding adsorption capacities of Cu 2+ and Zn 2+ are listed in Table 2.  www.nature.com/scientificreports/ The response model was obtained by fitting the experimental results with quadratic multiple regression equation. The response surface methodology (RSM) was used to construct the response relationship between Y q(Cu) , Y q(Zn) (Y q(Cu) and Y q(Zn) representing the adsorption capacity of biochar for Cu 2+ and Zn 2+ respectively) and preparation conditions (X 1 , X 2 and X 3 represent pyrolysis time, pyrolysis temperature and heating rate, respectively), and the interaction between preparation conditions and their influence on adsorption capacity were studied. The response models are shown in Eqs. (6) and (7), and in the model, the positive sign (+) represent the synergistic effect and the negative sign (−) indicates antagonistic effect 19 .
Analysis of variance (ANOVA) was performed on the experimental results of BBD, and the analysis results were shown in Table 3.
From Table 3, it can be found that the F-value for individual term of X 1 , X 2 and X 3 are 0.66, 28.97 and 2.43 for Cu-BC, indicating that the influence of preparation conditions on Cu 2+ adsorption performance of biochar was as follows: pyrolysis temperature (X 1 ) > heating rate (X 2 ) > pyrolysis time (X 3 ), and 6.22, 64.73 and 2.91 for Zn-BC, and the results showed that the influence of preparation conditions on Zn 2+ adsorption performance of biochar was as follows: pyrolysis temperature (X 1 ) > pyrolysis time (X 2 ) > heating rate (X 3 ). As can be seen from the above results, among the three preparation parameters, pyrolysis temperature has the greatest influence on the adsorption efficiency of Cu-BC and Zn-BC. The conclusion is consistent with many research, the pyrolysis temperature is the most effective pyrolysis factor to determine the adsorption capacity of biochar [20][21][22] . However, the sequence of the effects of pyrolysis time and heating rate on the two biochars was different. According to the F-value, the order of interaction items influence on adsorption performance of Cu-BC were 23 . The p-value < 0.05 indicates the significance of terms, less than 0.01, shows that a model term was considered as extremely significant [24][25][26] . In this study, the p-value of models for adsorption capacities of Cu 2+ and Zn 2+ were 0.0002 and < 0.0001, and all less than 0.01, indicated that two models were extremely significant. In this case, according to the p-value, the significant model terms for Cu 2+ removal are X 2 , X 1 X 2 , X 1 2 , X 2 2 , X 3 2 , and for Zn 2+ removal are X 1 , X 2 , X 1 2 , X 2 2 , X 3 2 . The coefficient R 2 and R adj 2 are used to further verify the validity of the model. In this study, the R 2 values for Eqs. (7) and (8) were 0.9662 and 0.9786, respectively. It indicates that about 96.62% and 97.86% of regression models were ascribed to the preparation conditions researched. In addition, the R adj 2 for Cu 2+ and Zn 2+ were 0.9229 and 0.9512, and the R 2 and R adj 2 of BC-Cu and BC-Zn were all above 0.92, indicating that the model has good accuracy. The Adeq. Precision (Adequate precision) is used to indicate the signal to noise ratio (SNR). The SNR > 4 is indicated that the signal is adequate, and the model can be used to guide the design space 27 . The SNR of Cu-BC and Zn-BC are 12.8699 and 16.6178, respectively, all much higher than 4, combined with the values of R 2 and R adj 2 , it showed that the two quadratic regression models are in good agreement with the experimental results.   Fig. 2. As can be seen from the figure, the predicted results of the model are close to the experimental values, indicating that the correlation between the preparation conditions of target biochars and target heavy metal ions adsorption capacity.

Effect of preparation conditions on adsorption capacities of target heavy metal ions by target biochars.
According to quadratic multiple regression model Eqs. (6) and (7), the response surface threedimensional diagram of Cu 2+ and Zn 2+ adsorption capacity of biochar can be obtained by the interaction of pyrolysis time (X 1 ), pyrolysis temperature (X 2 ) and heating rate (X 3 ). In Fig. 3 (Cu) and Fig. 4 (Zn), the interaction of the pyrolysis time, pyrolysis temperature, and heating rate on the adsorption capacity of target heavy metal ions can be observed.
For Cu-BC, as shown in Fig. 3a,b, the increment of adsorption capacity of Cu-BC for Cu 2+ occurred from 2.0 until 3.5 h of the pyrolysis time, and it reached the optimum time about 3.1 h. In terms of pyrolysis temperature, the adsorption capacity of Cu 2+ increased from 200 to 450 °C. It can be obtained in Fig. 3a,c, the optimum pyrolysis temperature at about 400-500 °C, the adsorption capacity of Cu 2+ by Cu-BC was increased, and the best pyrolysis temperature about 430 °C for Cu-BC. Meanwhile, Fig. 3b,c showed the heating rate for preparing Cu-BC, and the optimum heating rate is 19.5 °C/min for Cu-BC.
In Fig. 4a-c, the optimum preparation conditions for Zn-BC were shown. As for pyrolysis time, the adsorption capacity of Zn 2+ increased gradually from 1.5 to 2.5 h, and the pyrolysis time was about 2.2 h when the maximum adsorption capacity was reached. The optimum pyrolysis temperature for Zn-BC was displayed in Fig. 4a,c, similar to the optimum pyrolysis temperature of Cu-BC, the optimum pyrolysis temperature of Zn-BC also occurs in the range of 400-500 °C, and the highest adsorption capacity was achieved about 420 °C. The heating rate for Zn-BC was shown in Fig. 4b,c, different from the Cu-BC, the optimum heating rate interval was 14-16 °C/min, and the best optimum heating rate for Zn-BC is 15.8 °C/min.

Optimization of preparation conditions and verification of results.
Through the optimization analysis of the experimental results by Design-Expert 10.0 software, the optimal combination of preparation conditions of two target biochars were obtained. For Cu-BC, the optimized parameters were pyrolysis time of 3.09 h, pyrolysis temperature of 425 °C and heating rate 19.65 °C/min, and as for Zn-BC, the optimized prepa-   Table 4 and the validation results are displayed in Table 5. It can be clearly seen that the predicted results (18.12 and 19.63 mg/g for Cu-BC and Zn-BC, respectively) are in good agreement with the actual experimental results (18.92 and 19.98 mg/g for Cu-BC and Zn-BC, respectively), and the error is small. The factor named "desirability" from software were used to showed the best optimum parameters for output, the value of the "desirability" is closer to 1, the best conditions can be achieved 19,27,28 . In Table 5, the desirability values for Cu-BC and Zn-BC are 0.935 and 0.999, respectively, were closer to 1, and indicated that the best parameter conditions can be achieved. Therefore, the software has given a good prediction for pyrolysis conditions for M-BCs (M means heavy metal ions), and the target biochar prepared under the optimized preparation conditions has a strong adsorption capacity for the target heavy metal ion, which is technically and economically feasible.

Target heavy metal ions adsorption kinetics of target biochars. The kinetics curves for Cu 2+
removal by BC and Cu-BC, and Zn 2+ removal by BC and Zn-BC are displayed in Fig. 5a,b (BC in here is the unoptimized biochar which prepared at 2 h, 400 °C and 20 °C/min). As can be seen from the Fig. 5, the adsorp-    [29][30][31][32] . This indicated that in the rapid adsorption stage, the adsorption of heavy metal ions by biochar mainly occurs on the surface of biochar, with physical adsorption predominating. After 1 h, the adsorption rate begins to slow down, and chemical adsorption controls the adsorption process until the equilibrium is reached, and in this stage, the heavy metal ions diffuse into the micropores of the biochar and bind to the inner surface of the biochar to reach the adsorption endpoint. The BC and M-BCs showed that a similar trend displayed in kinetics, at the beginning of adsorption, the adsorption rate is fast, and about 80% of the adsorption removal occurred in first 1 h. The adsorption of BC on Cu 2+ and Zn 2+ reached the adsorption moderation point at about 50 min, and he adsorption equilibrium of BC for the two heavy metal ions was reached at about 120 min, while the target biochar (M-BCs) had a high adsorption rate in the first 30 min, and continued to maintain a certain adsorption level in the following time, and reached the equilibrium at about 4 h. However, in the adsorption kinetics experiment, the concentration of heavy metal ions selected was higher than that of BC due to the strong adsorption capacity of targeted biochar (M-BCs). As a result, the saturated adsorption capacity of the two biochars for heavy metal ions was greatly different, and the discussion on adsorption performance was more inclined to adsorption isotherm.
Based on the experimental results, the kinetic models were fitted, and the constants and coefficients obtained are listed in Table 6. It can be seen from Table 6, after fitting convergence, the correlation is extremely significant (R 2 > 0.9), indicating that the adsorption of heavy metal ions by BC and M-BCs included physical adsorption and chemical adsorption. However, compared to the pseudo-first-order kinetic, as shown in Table 6, Cu 2+ adsorption by BC and Cu-BC, Zn 2+ adsorption by BC and Zn-BC fitted better to the pseudo-second-order kinetic (R second 2 > R first 2 ), which explained that the pseudo-second-order model preferably can well define the adsorption process of biochar for heavy metal ions. The pseudo-second-order model shows that the adsorption is mainly controlled by chemical action, not by material transport steps, and the pseudo-second-order model of reaction adsorption rate is mainly based on surface chemical reaction control, such as surface complexation and precipitation 32 .

Target heavy metal ions adsorption equilibrium isotherms of target biochars.
In order to further prove that the M-BCs had better adsorption performance than that of the BC, the Langmuir and Freundlich isotherms models were used to fit the adsorption equilibrium results of BC and M-BCs for Cu 2+ and Zn 2+ . The fitting curves and parameters of the two models are displayed in Fig. 6 and Table 7, respectively.
As can be seen from the Fig. 6a,b, the adsorption effect of Cu-BC on Cu 2+ , Zn-BC on Zn 2+ are significantly better than that of BC. The fitting results of Langmuir model are similar to the measured values, indicating that the fitting reliability of the adsorption results is high. According to the fitting parameters in Table 7, the correlation coefficient R 2 , except for adsorption of Cu 2+ by BC, the Langmuir model of Cu-BC, BC (Zn 2+ ) and Zn-BC had a better fitting than the Freundlich model, which illustrated that a monolayer type adsorption onto  www.nature.com/scientificreports/ the homogenous sites of biochars 16,33 . The maximum adsorption capacities (Q max ) were fitted by the Langmuir model of BC (Cu 2+ ), Cu-BC, BC (Zn 2+ ) and Zn-BC were 177.66, 210.56, 146.14 and 223.32 mg/g. It was further proved that the adsorption performance of target biochars for targeted heavy metal ions was obviously stronger than that of the unoptimized biochar for heavy metal ions.

Characterizations of biochars.
The Zeta potential of BC, Cu-BC and Zn-BC are shown in Fig. 7. As shown in Fig. 7, the pH pzc of the BC, Cu-BC and Zn-BC were 2.44, 2.70 and 2.62. When the pH value > 2.70, the BC, Cu-BC and Zn-BC all have negative charges on the surface 23 , and there will be strong electrostatic attraction between them and the positive charges on the surface of heavy metal ions 34 . At the room temperature of 25 °C, the solubility product constant (K sp ) of Cu(OH) 2 and Zn(OH) 2 are 5.0 × 10 -20 and 7.1 × 10 -18 , respectively. According to the K sp of metal hydroxide, the pH of metal hydroxide precipitation at different concentrations can be calculated. The calculation formula is shown in Eq. (8): where, K sp and K w are the solubility product constant of metal hyoxide and ion product constant of water, respectively;α M n + is the concentration of metal ions (mol/L), Where M represents the metal ion and n represents the valence of the metal ion. When the reaction temperature and the type of heavy metals are determined, according to Eq. (8), the pH that affects the precipitation of heavy metal hydroxide is only the concentration of heavy metal ions, and the higher the concentration, the lower the pH. In this study, the concentrations of Cu 2+ and Zn 2+ solutions used in the experiment all were 50 mg/L, according to the Eq. (8), it can be calculated that the pH of Cu(OH) 2 and Zn(OH) 2 precipitation are 5.90 and 6.98, respectively. The pH of BC, Cu-BC and Zn-BC are 9.33, 9.73 and 9.55, which are all higher than the pH of the precipitation of the two metal hydroxides. Therefore, precipitation is one of the mechanisms by which biochar adsorb heavy metal ions.
The SEM images of Cu-BC and Cu-BC + Cu (the Cu-BC after adsorption of Cu 2+ ) were shown in Fig. 8a,b. Figure 8c,d displayed the SEM images of Zn-BC and Zn-BC + Zn (the Zn-BC after adsorption of Zn 2+ ). It can be seen from the SEM images that surface of two biochars were heterogeneous, cracked, and there are a large number of fragments. These accumulated fragments make the biochar form pores, and these pores enable heavy metal ions to enter the biochar. In order to further determine the elemental composition of biochar after adsorption, the EDX analysis of biochar was carried out.   35 . So, these metals can be exchanged by heavy metal ions in solution by cation exchange or coprecipitation of surface complexes during the adsorption process, and the heavy metal ions were transferred from aqueous solution to biochar 36,37 .
The analysis of specific surface area pore size distribution of BC, Cu-BC and Zn-BC are presented in Table 8 and Fig. 10a, the pore diameter of BC, Cu-BC and Zn-BC are 14.505, 21.686 and 8.357 nm, which are belong to mesoporous materials and are suitable for adsorption materials 38 . But the specific surface area of Cu-BC and Zn-BC are all smaller than the BC, and this is inconsistent with the view that the specific surface area is large and the adsorption capacity is enhanced, and this also further indicates that the performance of biochar obtained by characterization means does not necessarily explain the adsorption performance of biochar. Therefore, the characterization of biochar can only be used as an auxiliary means to explain the adsorption mechanism of biochar to adsorbents, but it is insufficient to measure the adsorption performance.
Nitrogen adsorption-desorption curves of the three adsorption materials are shown in Fig. 10b, the curves of BC, Cu-BC and Zn-BC are conformed to the IV isotherm adsorption, and type IV isotherms are derived from mesoporous adsorbent materials, with multilayer adsorption and capillary condensation phenomenon, and the hysteresis loop is H3 model, showed that the adsorption channel is relatively narrow, and the adsorption aperture for sheet pile holes 39 .
The functional groups of biochar before and after adsorption were analyzed by FTIR, and the FTIR spectra of BC, Cu-BC and Cu-BC + Cu, Zn-BC and Zn-BC + Zn are illustrated in Fig. 11a-c. It can be seen from the Fig. 11a, the peak pattern of BC, Cu-BC and Zn-BC are almost similar, there are only some difference were noticed in the intensity and wavenumbers, which means there are only difference in the surface chemistry, and the spectra were consistent with the water hyacinth biochar which prepared by Hashem et al. 40 . Peaks around 3450/cm was observed for BC, Cu-BC and Zn-BC, attributed to the stretching of -OH, this peak exists in most biochars are prepared from water hyacinth. The transmittance of BC near 2368/cm was highest in bichars, and the intensity in the Cu-BC and Zn-BC decreased. Peaks around 1648/cm for biochars were indicated to the aromatic C=C. A peak at 1084/cm is attributed to the stretching of C-O-C stretch associated with -OH bending of cellulose, hemicellulose, and lignin. 874/cm is due to aromatic C-H out of plane bends [41][42][43] . The FTIR patterns after adsorption of Cu 2+ and Zn 2+ are displayed in Fig. 11b,c, the width and location of these peaks changed, especially at 3400 and 1500/cm, it indicates that O-H groups plays an important role in the adsorption of heavy metal ions by biochar. Compared with before adsorption, the spectra of Cu-BC + Cu and Zn-BC + Zn showed obvious peaks at 1637 and 1639/cm of the carbonyl functional group (C=O), indicating the presence of metallic carbonyl bonds 29 .      Figure 12b,c represented the XRD pattern of Cu-BC and Cu-BC + Cu, Zn-BC and Zn-BC + Zn, respectively. From Fig. 12a, there are more peaks were observed in the XRD results of Cu-BC and Zn-BC than in BC, which revealed that the optimized biochar contained more mineral components than the common biochar did 44 . As www.nature.com/scientificreports/ can be seen from Fig. 12b,c, after adsorption of Cu 2+ and Zn 2+ , four major peaks of Cu-BC and Zn-BC all disappeared. For Cu-BC, after adsorption of Cu 2+ , the diffraction peak at 2θ = 29.48° shifts to 29.26°, and the intensity decreases, which is also the case for Zn-BC (2θ = 29.64° of Zn-BC and 29.34° of Zn-BC + Zn). These results indicate that when biochar reacts with Cu or Zn, Cu or Zn compounds are produced. Combined with EDX and FTIR analysis, Cu 2+ or Zn 2+ reacts with OH − and CO 3 2− in biochar to generate precipitate, and ion exchange reactions can also occur with cations K, Mg and Ca in biochar. It was successfully exchanged out of the water by biochar and to achieve the purpose of removal.
Through the characterization of biochar before and after adsorption, combined with the corresponding chemical basis, the adsorption mechanism of biochar for heavy metal ions mainly includes precipitation reaction, surface physical adsorption, cation exchange, electrostatic adsorption and surface complexation, etc.

Conclusions
In this study, the effects of preparation conditions: pyrolysis time, pyrolysis temperature and heating rate on the performance of target biochar for target heavy metal ions adsorption was studied. The RSM coupled with BBD was applied to optimize the preparation parameters of target biochars. With Cu 2+ and Zn 2+ as the target pollutants, the production conditions of biochar adsorbed on the two heavy metal ions were optimized. The optimization results reported that two different optimal combination conditions for biochar with adsorption of Cu 2+ and Zn 2+ . For Cu-adsorbed biochar, the optimal production conditions are: pyrolysis time of 3.09 h, pyrolysis temperature of 425 °C, and heating rate of 19.65 °C/min. While for Zn-adsorbed biochar, the optimal production conditions are: pyrolysis time of 2.19 h, pyrolysis temperature of 422 °C, and heating rate of 15.88 °C/min. The adsorption performance of Cu-BC and Zn-BC to their target pollutants is better than that of the biochar without optimized production conditions. The Langmuir isothermal model fitting the maximum theoretical adsorption capacities of BC and Cu-BC on Cu 2+ are 177.66 and 210.56 mg/g, respectively, and the maximum theoretical adsorption capacities of BC and Zn-BC on Zn 2+ are 146.14 and 223.32 mg/g, respectively. As general conclusion, the optimized target biochars Cu-BC and Zn-BC are effective adsorbents, which can be used for the removal of target heavy metals. This study not only prepared targeted biochar for target pollutant, but also reduced the number of experiments and the number of samples characterized to a certain extent, which has a great effect on saving raw materials, time and production costs.